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Abstract: Land use/cover change (LUCC) is becoming more and more frequent and extensive as a result 
of human activities, and is expected to have a major impact on human welfare by altering ecosystem 
service value (ESV). In this study, we utilized remote sensing images and statistical data to explore the 
spatial-temporal changes of land use/cover types and ESV in the northern slope economic belt of the 
Tianshan Mountains in Xinjiang Uygur Autonomous Region, China from 1975 to 2018. During the study 
period, LUCC in the study region varied significantly. Except grassland and unused land, all the other land 
use/cover types (cultivated land, forestland, waterbody, and construction land) increased in areas. From 
1975 to 2018, the spatial-temporal variations in ESV were also pronounced. The total ESV decreased by 
4.00108 CNY, which was primarily due to the reductions in the areas of grassland and unused land. 
Waterbody had a much higher ESV than the other land use/cover types. Ultimately, understanding the 
impact of LUCC on ESV and the interactions among ESV of different land use/cover types will help 
improve existing land use policies and provide scientific basis for developing new conservation strategies 
for ecologically fragile areas. 
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1 Introduction 


Land use/cover change (LUCC) is a significant component of global environmental change. 
This is because LUCC reflects interactions between human activities and natural resources. 
LUCC is driven by multiple factors and has both direct and indirect impacts on environment, as 
well as on global sustainable development. This makes LUCC being one of the most 
decisive factors affecting global ecosystems (Xiao et al., 2016; Sun et al., 2019). Range and 
strength of LUCC ina given region often directly influence primary ecological processes 
and thus ecosystem functioning (Long et al., 2014; Arowolo and Deng, 2018; Zhang et al., 
2018; Cui et al., 2019). The rapid growth of population and expansion of human activities may 
hinder the provision and the maintenance of vital ecosystem services (Paz-Kagan et al., 2014; 
Li et al., 2018; Wang et al., 


“Corresponding authors: MA Yonggang (E-mail: mayg@xju.edu.cn); GONG Lu (E-mail: gonglu721@163.com) 

Received 2021-03-18; revised 2021-07-02; accepted 2021-07-06 

© Xinjiang Institute of Ecology and Geography, Chinese Academy of Sciences, Science Press and Springer-Verlag GmbH Germany, 
part of Springer Nature 2021 


http://jal.xjegi.com; www.springer.com/40333 


SUN Chen et al.: Response of ecosystem service value to land use/cover change in the... 


2019). The ecosystem service value (ESV) concept, which puts a monetary value on the 
ecosystem services provided by natural areas, was first proposed by Costanza et al. (1997). This 
approach provides an assessment on the monetary value of ecological processes which humans 
depend on and are closely related to human well-being. Ecosystem services include supplying 
services, regulating services, supporting services, and cultural services (Xiao et al., 2016; Chen et 
al., 2019). Numerous studies have shown that LUCC significantly influences ESV, as well as the 
interactions of ESV among different ecosystem service types and different ecosystem types (Cord 
et al., 2017; Huang et al., 2021; Liu et al., 2021; Pan et al., 2021). Therefore, a better 
understanding of how LUCC affects ESV and its interactions with time is essential for the 
development of new ecosystem management policies (Huang et al., 2019). 

The ESV approach has become an essential tool for evaluating ecosystem services across 
ecosystems (Costanza et al., 2014; Yuan et al., 2019). Since 1997, many studies have been 
conducted on the impact of LUCC on ESV from local to global scales (Nahuelhual et al., 2014; 
Wang et al., 2014). At present, ESV has been quantified for different ecosystem types, such as 
forests (Xie et al., 2010; Li et al., 2017), grasslands (Homolova et al., 2014), wetlands (Wang et 
al., 2006; Gunderson et al., 2016), farmlands (Kroeger and Casey, 2007), and marine 
environments (Barbier et al., 2011). As a result, many ESV assessment methods have been 
developed, including cost-benefit analysis (de Groot et al., 2012), willingness to pay assessment 
(Rawlins and Westby, 2013), value transfer (Xie et al., 2015), public participation in geographic 
information system (GIS) method (Brown, 2013), selection experiment (Jobstvogt et al., 2014), 
and conditional valuation (Brander et al., 2012). In China, ESV studies have been focused on 
ecosystem service regimentation, value calculation, quantitative assessment, and spatial identity 
(Ouyang et al., 2016; Xiao et al., 2016). For example, Cao et al. (2018) assessed the gains and 
losses of ESV, which triggered by the rapid urbanization, in the coastal areas of Zhejiang Province 
from the LUCC perspective; Wang et al. (2018) identified the elements promoting LUCC and 
quantified the relationship between ESV and elevation in the Hengduan Mountains. In the past 
few years, with the development of remote sensing, GIS, and global positioning system (GPS) 
technologies, ESV research was increasingly focused on assessing individual ESV for multiple 
ecosystem types (Ding et al., 2020; Liu et al., 2021). This strategy has promoted the quantitative 
evaluation of ESV and also has encouraged the analysis of spatial and temporal patterns for ESV. 
However, few studies considered how spatial heterogeneity in LUCC affects ESV, nor offered a 
visual representation of these spatial changes (Yuan et al., 2019). 

China has experienced significant economic development since 1978, especially in the eastern 
developed regions such as the Yangtze River Delta and the Pearl River Delta. As a result of 
economic development, LUCC has accelerated in the whole country (Chen et al., 2019; Hu et al., 
2019; Liu et al., 2019; Zhang et al., 2019). Rapid economic development and urban expansion 
have caused significant changes in land use/cover, which have in turn resulted in changes of 
ecosystems near cities and changes in ESV. In particular, urban growth in arid regions occurred at 
the expense of reduced farmland and ecological stability under threatened. Despite these 
challenges, detailed and in-depth studies on the quantitative effects of how LUCC influences the 
provision of ecosystem services in arid regions remain lacking (Tang et al., 2017). Therefore, new 
research on relationship between LUCC and ecosystem services in arid regions is vital for 
providing theoretical basis for land management policy making (Huang et al., 2019). 

The northern slope economic belt of the Tianshan Mountains is situated in the middle of the 
Xinjiang Uygur Autonomous Region in Northwest China, which is part of the core area of the 
Silk Road Economic Belt. The northern slope economic belt of the Tianshan Mountains is also the 
primary location of agriculture and human settlement in Xinjiang (Sun et al., 2018; Wang et al., 
2018). Over the past 44 a, rapid urbanization has led to a dramatic shift in land use patterns in this 
region, affecting the provision of ecosystem services (Li et al., 2019). 

Here, we quantified ESV of the northern slope economic belt of the Tianshan Mountains from 
1975 to 2018, and used land use/cover transfer matrix to evaluate the temporal and spatial 
patterns of LUCC. The purposes of this study included: (1) quantifying the temporal and spatial 
dynamics of LUCC within the study region from 1975 to 2018; (2) assessing the changes of ESV 
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and identifying the related shifts of land use/cover; and (3) exploring the interactions between 
LUCC and ecosystem services. 


2 Materials and methods 


2.1 Study area 


The northern slope economic belt of the Tianshan Mountains (40°43'—47°15'N, 79°53'—96°23’E) 
is situated in the hinterland of the Eurasian continent, along the southern edge of the 
Gurbantunggut Desert, China. The region covers an area of 9.54x10* km’, accounting for 5.7% of 
the total area of Xinjiang. The elevation of this region is higher in the southeast and lower in the 
northwest, ranging from 400 to 600 m. The climate varies from temperate arid to semi-arid, with 
the annual mean temperature of 6.9°C, mean annual precipitation of 220 mm, and mean annual 
evaporation of 1817 mm. This area experiences a range of adverse weather conditions, including 
strong winds and sand storms. According to the China Statistical Yearbook (National Bureau of 
Statistics, 2020), the total human population in this region was 45.8x 10° in 2018, accounting for 
23.3% of the total population in Xinjiang. The gross domestic product (GDP) was about 
8.51x10'? CNY in 2018. With the implementation of the Belt and Road Initiative, the northern 
slope economic belt of the Tianshan Mountains becomes a critical development hub in Central 
Asia; it is also part of the core area of the Silk Road Economic Belt (Sun et al., 2018; Wang et al., 
2018). 


2.2 Data collection 


We used land use/cover data and socio-economic data during the period of 1975-2018 (Local 
Chronicles Compilation Committee of Xinjiang Uygur Autonomous Region, 1975-2018). Land 
use/cover data were obtained from remote sensing images, which were taken on 2 August 1975 
(Landsat MSS), 20 July 1990 (Landsat TM), 12 July 2005 (Landsat ETM+), and 15 July 2018 
(Landsat TM) (http://www.gscloud.cn/). Geometric correction of images was performed using the 
1:100,000 National Standard Terrain Database. Errors were controlled within 0.5 pixels. We 
classified all images by combining supervised classification with visual inspection and 
re-sampled them at a spatial resolution of 30 m. We used remote sensing images to identify land 
use/cover types in the research area (General Administration of Quality Supervision et al., 2017). 
Land use/cover types included the following classifications: cultivated land, forestland, grassland, 
waterbody, construction land, and unused land (Zhang et al., 2019). We consulted local land use 
maps and conducted ground truth surveys to ensure the accuracy of the classifications. We 
selected 100 randomly generated sample points to conduct field surveys and made corrections to 
verify the actual land use distribution. In the end, the overall accuracy of each map was higher 
than 85%. The accuracy of remote sensing land use classification was estimated to be >85% using 
a subset of regional data from which real ground data can be obtained. 


2.3 Methods 


2.3.1 Dynamics of LUCC 


The K index quantitatively describes the dynamics of LUCC (Ye et al., 2018), which is calculated 
as follows: 


eee) 00%, (1) 
U T 


where K measures the land use dynamics during the study period (%); U, and Um respectively 
indicate the area of a certain land use/cover type at the beginning of the study period and at the 
end of the study period (km’); and T is the length of the study period (a). 

Land use transfer matrix can be used to capture the dynamics of LUCC for a given region over 
a specific period of time. It takes the form of a two-dimensional array populated by analyzing 
land use/cover images over time. Land use transfer matrix is ordinarily described as follows (Yan 
et al., 2019): 
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where Sj represents the percentage of the total land area converted from land use/cover type i to 
land use/cover type j during the study period (%); and S,, represents the percentage of area of 
land use/cover type n remained unchanged during the study period (%). 


2.3.2 ESV 


We estimated ESV using methods described by Costanza et al. (2014) and Xie et al. (2003), who 
had already described the ESV characteristics of China's ecosystems by consulting 200 Chinese 
ecologists. In order to quantity ESV, we introduced the concept of ESV equivalent value 
coefficient in this study to describe the money value of ESV for each ecosystem service function 
to each land use/cover type in unit time and area (i.e., per hectare per year). According to previous 
researches, we set ESV equivalent value coefficient as the annual output value of natural crop 
production per hectare, which is one-seventh of the actual annual output value of crop production 
(Xie et al., 2005). The average actual annual output value of crop production between 1975 and 
2018 was 3029.15 CNY/(hm~.a). 

Table 1 shows ESV equivalent value coefficients for different ecosystem service functions 
under different land use/cover types in the northern slope economic belt of the Tianshan 
Mountains (Xie et al., 2005). Based on Table 1, ESV for each ecosystem service function (ESV;) 
was calculated by multiplying the area of each land use/cover type with ESV equivalent value 
coefficient for the corresponding land use/cover type, then summing the outcomes of all land 
use/cover types (Eq. 3). The total ESV was then calculated as the sum of ESV; across all 
ecosystem service functions (Eq. 4). 


ESV, = $ (a x VC, ), (3) 
i=l 


ESV = X'$ (4 xC), (4) 


i=l j=l 
where ESV; is the estimated ESV of ecosystem service function j in the study area (CNY); A; is 
the area of land use/cover type i (hmô; VC; is the ESV equivalent value coefficient, i.e., the 
equivalent value coefficient of ecosystem service function j to the corresponding land use/cover 
type i (CNY/(hm’-a)); ESV is the total ESV; for the study area (CNY); n is the number of land 
use/cover type; and m is the number of ecosystem service type. 
2.3.3 Cross-sensitivity coefficient 
We tested the sensitivity of the variations of ESV equivalent value coefficient as follows: 
(ESV; - ESV, )/ESV, 

(VC,, ~ VC, )NCu 

We calculated the percentage change of ESV within a range of +50% of ESV equivalent value 
coefficient. CSS, is the cross-sensitivity coefficient (CSS) of an ecosystem service function to the 
corresponding land use/cover type k in time period t. ESV, and ESV; are the total ESVs before 
and after the variation of ESV equivalent value coefficient (CNY), respectively. ESV, is the total 
ESV in time period t. VC, and VC, are the ESV equivalent value coefficients before and after 
the +50% change (CNY/(hm/.a)), respectively. The closer the value of CSS, approaches to zero, 
the less sensitive the estimated ESV is to the uncertainty of ESV equivalent value coefficient. 


(5) 


CSS, = 


3 Results 
3.1 Dynamics of LUCC 


Land use/cover types and their dynamics over time in the northern slope economic belt of the 
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Table 1 Ecosystem service value (ESV) equivalent value coefficient for each ecosystem service function to the 
corresponding land use/cover type in the northern slope economic belt of the Tianshan Mountains 


ESV equivalent value coefficient (CNY/(hm’sa)) 


Ecosystem Ecosystem iG J - F 
service type service function ea ae Forestland Grassland Waterbody CoE o Unuse 
land land land 
. Grain production 1756.90 175.69 527.07 175.69 0.00 17.57 
Supplying R tenal 
services EEE 175.69 4567.96 87.85 17.57 0.00 0.00 
production 
Gas regulation 878.45 6149.18 1405.53 0.00 0.00 0.00 
, cana 1563.65 1405.53 1581.22 808.18 0.00 0.00 
Regulating regulation 
eee Hydrological 1054.15 5622.11 1405.53 17,300.00 -4232.60 52.71 
regulation 
Waste treatment 2881.33 2301.55 2301.55 0.00 —18,548.09 17.57 
. Soil conservation 2565.09 6851.95 3425.97 17.57 0.00 35.14 
Supporting Biodiversit 
services ny 1247.41 5727.52 1915.03 4374.70 0.00 597.35 
conservation 
ane Asstheng 17.57 2248.84 70.28 7624.99 0.00 17.57 
services landscape 


Note: Negative value means that the ecosystem service function negatively affects land use/cover type. 


Tianshan Mountains are presented in Figure 1. Grassland, cultivated land, and unused land were 
the primary land use/cover types during 1975-2018. Indeed, these three land use/cover types 
made up >90% area of the northern slope economic belt of the Tianshan Mountains. From 1975 to 
1990, the areas of cultivated land, forestland, waterbody, and construction land were increased. In 
contrast, the proportion of grassland and unused land were decreased during 1975—1990. The area 
of construction land during 1975-1990 increased from 939.88 to 1255.53 km’, with an average 
increasing rate of 2.23% (Fig. 1). From 1990 to 2005, the areas of cultivated land, waterbody, and 
construction land increased. From 2005 to 2018, the areas of forestland, grassland, waterbody, 
and unused land decreased, while cultivated land and construction land expanded. Of all the land 
use/cover types, the area of construction land continually increased during the last 44 a, showed 
the most positive change (Fig. 2). However, the initial area of construction land was small, so the 
total area of construction land increased by mere 1328.34 km’. Cultivated land also became more 
abundant, and waterbody resources were added to meet irrigation demands. Although grassland 
was the most common land use/cover type in the study area, it decreased significantly during the 
study period. The area of unused land changed little from 1975 to 2018. 


(a) 1975 N (b) 1990 


(c) 2005 (d) 2018 


Legend 
B Cultivated land MM Forestland {| Grassland 0 325 650 km 
E Waterbody WH Construction land [_] Unused land 


Fig. 1 Distributions of different land use/cover types in the northern slope economic belt of the Tianshan 
Mountains in 1975 (a), 1990 (b), 2005 (c), and 2018 (d) 
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D Cultivated Forestland Grassland Waterbody Construction Unused 


land land land 


Fig. 2 Dynamics of land use/cover change (LUCC) in the northern slope economic belt of the Tianshan 
Mountains from 1975 to 2018 


Land use transition matrix showed that grassland and unused land were converted to other land 
uses during the study period (Table 2). For example, about 3047.78 km? of other land use/cover 
types were transformed to forestland. Specifically, from 1975 to 2018, the expansion of forestland 
was primarily converted from grassland, with the converted area of 1775.49 km?; soil and water 
conservation projects likely explained the progress of grassland conversion. A further 441.67 km? 
of the expansion was from unused land. 


Table 2 Land use transfer matrix for the northern slope economic belt of the Tianshan Mountains from 1975 to 


2018 
Period Land use/ Cultivated Forestland Grassland Waterbody Construction Unused Total 
cover type land (km’) (km) (km’) (km?) land (km) land (km’) (km?) 
Eo 13,100.00 204.36 2228.15 41.75 188.98 461.45 16,224.69 
Forestland 109.87 5324.26 202.16 6.59 13.18 39.55 5695.61 
1975- Grassland 2149.04 503.20 60,300.00 101.08 210.95 1052.55 643,16.82 
1990 Waterbody 32.96 6.59 59.33 4698.01 4.39 94.49 4895.77 
ee 127.45 2.20 32.96 0.00 711.95 8.79 883.35 
Unused land 257.09 50.54 1995.23 175.79 61.53 53,900.00 56,440.18 
a 5985.68 1758 2032.58 153.82 320.82 219.73 8730.21 
Forestland 10.99 4983.67 81.30 2.18 17.58 2.38 5098.10 
1990- Grassland 28.57 68.12 65,100.00 90.09 238.73 318.62 65,844.13 
2005 Waterbody 0.00 2.19 10.99 4384.39 13.18 10.98 4421.73 
ag 2.20 2.19 21.97 4.39 1153.63 6.59 1190.97 
Unused land 10.99 188.98 628.45 105.47 68.12 54,500.00 55,502.01 
oe 17,300.00 0.00 202.16 13.18 305.44 50.54 17,871.32 
Forestland 92.29 5915.36 6.59 4.39 19.78 0.00 6038.41 
2005- Grassland 1779.88 13.18 60,500.00 35.16 261.49 35.16 62,624.87 
2018 Waterbody 4.39 0.00 21.97 5216.59 3.26 17.58 5263.79 
Gie ai 17.58 0.00 2.21 8.79 1481.04 0.00 1509.62 
Unused land 1070.13 2.20 85.70 32.96 136.24 538.00 1865.23 
pe 36,385.68 221.94 4462.89 208.75 815.24 731.72 42,826.22 
Forestland 213.15 16,223.29 290.05 13.16 50.54 41.93 16,832.12 
1975- Grassland 3957.49 584.50 1859.00 226.33 711.17 1406.33 8744.82 
2018 Waterbody 37.35 8.78 92.29 14,298.99 20.83 123.05 14,581.29 
oa 147.23 4.39 57.14 13.18 3346.62 15.38 3583.94 
Unused land 1388.21 241.72 2709.38 314.22 265.89 1622.00 6541.42 
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Waterbody was mostly converted to grassland and unused land. Meanwhile, 1328.57 km’ of 
construction land was converted to other land use/cover types, and most of it was converted to 
cultivated land. Overall, the area of construction land increased by 3.29% during the study period 
(Fig. 2), primarily owing to the conversion of cultivated land and unused land. The next largest 
gain was cultivated land (0.58%), followed by waterbody (0.20%) and forestland (0.12%). These 
increases were compensated by the reductions of grassland and unused land. Hence, given these 
trends (i.e., areas increased in construction land and cultivated land), there might be growing 
demands on soil water for vegetation restoration in the study area. 


3.2 Changes in ESV 


3.2.1 Variations in ESV among different land use/cover types 


As shown in Figure 3 and Table 3, the total ESV of the study region decreased over the study 
period. Overall, grassland was the major contributor to ESV among the six land use/cover types. 
Indeed, it accounted for 50.52%—53.33% of the total ESV. ESV of grassland increased slightly in 
2005, but then fell by 2.23% during 2005-2018. This decline was caused by the large scale 
conversion of grassland to cultivated land and construction land. Waterbody made the second 
highest contribution to the total ESV (19.81%—21.85%), even though its area only accounted for 
3.31%-3.57% of the study region. Over the entire study period, the proportion of waterbody to 
the total ESV gradually increased (from 19.81% in 1975 to 20.48% in 2018), while the 
contribution of grassland declined from 53.33% to 50.52%. As seen in Table 3, the mean annual 
rates of ESV change for cultivated land, forestland, and construction land were all greater than 
zero, indicating that ESV of these land use/cover types increased over time. This increase was 
mostly marked by construction land, at the rate of 3.29% (Fig. 2). In contrast, ESV of 
construction land was negative (Table 3). Mean annual rates of ESV change for unused land were 
all lower than zero (Table 4), indicating an overall reduction in ESV of this land use/cover type. 


9) mm 1975 
E 1990 
mm 2005 
lm 2018 


Mi i W W 


Cultivated Forestland Grassland Waterbody Construction Unused 
land land land 


ESV (x10!°CNY) 
or N LU Hh ADs w 


| 
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Fig. 3 Ecosystem service value (ESV) for different land use/cover types in the northern slope economic belt of 
the Tianshan Mountains in 1975, 1990, 2005, and 2018. Negative value means that the loss caused by 
depreciation of natural capital and environmental degradation. 


For waterbody, the largest increase in ESV was observed during 1975-1990, obtaining a gain 
of 2.70x10° CNY, which constituted the maximum mean annual rate of ESV change at 1.80% 
(Table 4). ESV of construction land decreased by 8.00x10° CNY during 1975-1990, declining by 
0.53%/a on average. From 1990 to 2005, the greatest gain in ESV was observed for cultivated 
land (2.20x 10° CNY), with a mean annual increasing rate of 1.47%. In contrast, ESV of grassland 
declined by 2.10x10° CNY, with —1.40% mean annual increasing rate on average. Cultivated land 
earned the largest growth in ESV from 2005 to 2018, with an increase of 2.7010’ CNY and a 
mean annual increasing rate of 2.08%. Waterbody decreased in ESV by 3.40x10° CNY, displaying 
a mean annual increasing rate of —2.62% in 2005-2018. Over the entire study period, the impact 
of LUCC on ESV differed among land use/cover types. Among them, cultivated land had the 
largest mean annual rate of ESV change. 
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Table 3 ESV for different land use/cover types in the northern slope economic belt of the Tianshan Mountains 
in 1975, 1990, 2005, and 2018 


Land use/cover type 


Year ESV Itivat tructi d 
Cultivated Forestland Grassland Waterbody Construction Unuse 


land land land 

ESV (x10'° CNY) 1.90 1.99 8.17 3.05 —0.21 0.42 

a Ga lc 12.41 12.98 53.33 19.91 -1.37 2.74 
ESV (x10 CNY) 1.97 2.16 8.22 3.32 —0.29 0.41 

1330 ka aaa 12.48 13.68 52.06 21.03 -1.84 2.59 
ESV (x10'° CNY) 2.19 2.13 8.01 3.47 -0.36 0.44 

one eRe en 13.79 13.41 50.44 21.85 -2.26 2.77 
ESV (x10'° CNY) 2.46 2.09 7.72 3.13 -0.52 0.40 

oe poe 16.10 13.68 50.52 20.48 -3.40 2.62 


Note: Negative value means that the loss caused by depreciation of natural capital and environmental degradation. 


Table 4 Dynamics of ESV for different land use/cover types in the northern slope economic belt of the Tianshan 
Mountains from 1975 to 2018 


Land use/cover type 


Time period ESV F i 
p Cultivated Forestland Grassland Waterbody Construction Unused 


land land land 
Charge of ESV 0.07 0.17 0.05 0.27 -0.08 -0.01 
(x10 CNY) 
1975-1990 Mean annual rate of 
ESV change (%) 0.47 1.13 0.33 1.80 ~0.53 -0.07 
Chanko ot EON 0.22 -0.03 -0.21 0.15 -0.07 0.03 
(x10 CNY) 
1990-2005 Mean annual rate of 
ESV change (%) 1.47 -0.20 -1.40 1.00 -0.47 0.20 
(tiang LESY 0.27 -0.04 -0.29 -0.34 -0.16 -0.04 
(x10! CNY) 
2005-2018 Mean annual rate of 
ESV change (%) 2.08 -0.31 -2.23 -2.62 -1.23 -0.31 
Charges of PSY 0.56 0.10 -0.45 0.08 -0.31 -0.02 
(x10 CNY) 
1975-2018 Mean annual rate of 
ESV change (°) 1.30 0.23 -1.04 0.19 -3.43 -0.05 


3.2.2 Comparisons of ESV among different ecosystem service functions 


ESV of each ecosystem service type was calculated from 1975 to 2018 using ESV equivalent 
value coefficient of each ecosystem service function to the corresponding land use/cover type 
(Fig. 4), and the total ESV decreased by 4.00x 10° CNY during this period. In 1975, the total ESV 
was 1.53x10'' CNY. At this time, regulating services had the highest value, accounting for 
52.03% of the total ESV, followed by supporting services (37.22%), supplying services (6.84%), 
and cultural services (3.91%). In 1990, the total ESV was 1.58x10'' CNY, an increase of 
1.47x10'° CNY compared with 1975. This increase was primarily attributable to the rise of 
regulating services (1.28x10'° CNY). Supplying, supporting, and cultural services also increased 
by about 5.00x10°, 1.00x10°, and 4.00x10° CNY, respectively. From 1990 to 2005, the total ESV 
increased by 9.00x10° CNY, primarily owing to an increase of supplying services. In contrast, 
regulating services decreased by 8.00x10* CNY, and cultural services also decreased by 1.00x10° 
CNY. The total ESV decreased from 1.58x10'' to 1.53x10'' CNY from 2005 to 2018; the most 
important factor was waste treatment. From 1975 to 2018, regulating services and supporting 
services decreased by 1.20x10° and 1.00x10* CNY in the study region, respectively; while 
cultural services and supplying services increased by 1.00x 10° and 8.00x 108 CNY, respectively. 
From 1975 to 2018, waste treatment accounted for the largest proportion of the total ESV in the 
study region; it valued 3.46x10'° CNY in 1975, increased to 3.57x10'? CNY in 1990, and 
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dropped to 3.32x10'° CNY in 2018 (Fig. 4). The next most significant contributor to the total 
ESV was soil conservation, which valued 3.02x10'° CNY in 1975 and 3.12x10'° CNY in 2005, 
and dropped to 3.02x10'° CNY in 2018. Raw material production, aesthetic landscape, and grain 
production made the lowest contributions to the total ESV, with individual contributions never 
exceeding 7.00x10° CNY. From 1975 to 2018, ESV of climatic regulation fluctuated, ultimately 
increased by 2.37% (Fig. 5). Grain production also increased ESV by 0.14% and this increase was 
positively correlated with the expansion of cultivated land. Waste treatment decreased ESV by 
0.33%, followed by hydrological regulation and biodiversity conservation. 
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Fig. 4 ESV for different ecosystem service functions in the northern slope economic belt of the Tianshan 
Mountains from 1975 to 2018 
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Fig. 5 Mean anneal rate of ESV change for different ecosystem service functions in the northern slope economic 
belt of the Tianshan Mountains from 1975 to 2018 


3.2.3 Spatial patterns of ESV 


We classified spatial patterns of ESV into five levels using natural breaks: lowest ecological 
environment risk (I), low ecological environment risk (II), moderate ecological environment risk 
(III), high ecological environment risk (IV), and highest ecological environment risk (V). High 
and highest ecological environment risk levels occurred primarily in construction land and 
waterbody, while low and lowest ecological environment risk levels were mainly observed in 
forestland, grassland, and cultivated land (Fig. 6). From 1975 to 2018, the area proportion of 
lowest and high ecological environment risk levels changed significantly (Table 5). In this period, 
the area with lowest ecological environment risk level decreased to 3.20x10* km’, and the area 
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with low ecological environment risk level reduced to 2.13x10* km’. Meanwhile, the area with 
high ecological environment risk level increased from 3.05x10* km? in 1975 to 3.14104 km? in 
1990, and rapidly expanded to 3.92x10* km? in 2005. Ecological environment risk in the study 
area was primarily affected by urbanization and LUCC. As a sensitive ecological region, the 
northern slope economic belt of the Tianshan Mountains is situated in an area with high 
ecological environment risk. In this study, the area with high ecological environment risk level 
was found to be increasing rapidly over time, while the area with lowest ecological environment 
risk level declined under the influences of economic development and urbanization. 
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Fig. 6 Spatial distributions of areas with different ecological environment risk levels in the northern slope 
economic belt of the Tianshan Mountains in 1975 (a), 1990 (b), 2005 (c), and 2018 (d). I, lowest ecological 
environment risk; II, low ecological environment risk; III, moderate ecological environment risk; IV, high 
ecological environment risk; V, highest ecological environment risk. 


Table 5 Variations of areas with different ecological environment risk levels in the northern slope economic belt 
of the Tianshan Mountains in 1975, 1990, 2005, and 2018 


Ecological 1975 1990 2005 2018 
environment Area Proportion Area Proportion Area Proportion Area Proportion 
tisk level (x 10* km’) (%) (x10* km’) (%) (x10* km’) (%) (x10* km’) (%) 
I 3.82 26.43 3.68 25.11 3.44 22.88 3.20 21.30 
II 2.81 19.41 2.88 19.65 2.14 14.20 2.13 14.20 
Ul 2.73 18.88 2.84 19.44 2.94 19.52 3.05 20.33 
IV 3.05 21.12 3.14 21.45 3.92 26.04 4.03 26.82 
Vv 2.05 14.16 2.10 14.35 2.61 17.36 2.61 17.35 


Note: I, lowest ecological environment risk; II, low ecological environment risk; III, moderate ecological environment risk; IV, high 


ecological environment risk; V, highest ecological environment risk. 


3.3 Sensitivity analysis 


We calculated CSS for ESV dynamics, driven by land use transformations in the study area, as 
shown in Figure 7. CSS represents all possible transitions among the six land use/cover types in 
the northern slope economic belt of the Tianshan Mountains. Since the net conversion area 
between two land use/cover types was symmetric (same for either type), we only showed one-way 
CSS for analysis. The most visible land use/cover type transition in the study area was from 
forestland to unused land, followed by the transformation from cultivated land, forestland, and 
grassland to waterbody. These changes therefore had the greatest impact on the total ESV. 
Although waterbody accounted for only about 3.50% of the total study area, its CSS was very 
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high. This is due to the fact that waterbody had a much higher ESV equivalent value coefficient 
than the other land use/cover types. 
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Fig. 7 Cross-sensitivity coefficient (CSS) of ESV for cultivated land (a), forestland (b), grassland (c) and 
waterbody (d) to represent the sensitivity of variations of ESV equivalent value coefficient for different land 
use/cover types in the northern slope economic belt of the Tianshan Mountains from 1975 to 2018 


4 Discussion 


The results indicate that the areas of the six land use/cover types in the northern slope economic 
belt of the Tianshan Mountains have changed significantly over the last 44 a (from 1975 to 2018). 
Land use patterns in the northern slope economic belt of the Tianshan Mountains were affected by 
changes of cultivated land area and infrastructure. With the intensification of urbanization and the 
further development of local economy, the amount of human infrastructure has increased 
dramatically (Huang et al., 2019). The northern slope economic belt of the Tianshan Mountains is 
the most densely populated area in Xinjiang. The expansion of urban areas in recent decades has 
coincided with the conversion of a large amount of cultivated land. Newly added infrastructure 
occurred mostly around major cities within the study area, and small- and medium-sized towns 
which surrounding these major cities were becoming more developed as well. The transformation 
of cultivated land to construction land was a consequence of urbanization and economic 
development, which has led to the occupation of previously cultivated land (Paz-Kagan et al., 
2014). Among all studied land use/cover types, the greatest change was observed for construction 
land. The two situations mentioned above made the expansion of these cities (Lausch et al., 
2015). In regard to the return of cultivated land to grassland, policy shifts encouraged the 
transition of cultivated land to forestland and grassland, therefore, a large quantity of cultivated 
land (scattered along the watercourse) was transformed to forestland and grassland from 1975 to 
2018. Compared with other regions in Xinjiang, the land area under cultivation with superior 
water resources has increased, likely due to the quality of agricultural land in the study region. 
The northern slope economic belt of the Tianshan Mountains is the most developed area of 
modern agriculture in Xinjiang. A growing regional population and high living standard have 
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increased the demand for agricultural products (Tang et al., 2017; Zhang et al., 2018). As a 
consequence, the amount of cultivated land throughout the study area was growing steadily over 
the last 44 a. 

LUCC can not only alter ecosystem structures and processes, but also influence regional 
ecosystems to provide products and services to human society (Paz-Kagan et al., 2014). ESV of the 
northern slope economic belt of the Tianshan Mountains decreased from 15.32 10'° CNY in 1975 
to 15.28x10'° CNY in 2018. This was the result of sacrificing ecological environment to regional 
economic development. The collective effect of natural factors and human activities leads to 
changes in ESV. Among these factors, LUCC is an essential driver of changes in ESV. ESV may be 
affected by the changes of land use structure, composition, pattern, and intensity. With the 
implementation of the Silk Road Economic Belt project, human activities, such as rapid economic 
development, population growth, accelerated urbanization, and land use policies, have had a 
significant effect on LUCC in the northern slope economic belt of the Tianshan Mountains (Lausch 
et al., 2015; Huang et al., 2019). 

The greatest proportion of ESV reduction was due to the conversion of grassland to construction 
land, which had a very low ESV equivalent value coefficient. In the last 44 a, the change of land 
use/cover in the northern slope economic belt of the Tianshan Mountains was significant. The 
dramatic increase in land used for human infrastructure was accompanied by a decrease in 
grassland. Large grassland area changed into cultivation land or construction land, subsequently 
resulting in a decrease in ESV. 

The method of ecosystem service evaluation proposed by Costanza et al. (2014) and the 
modifications proposed by Xie et al. (2003, 2005) have been shown to be reliable in many studies 
(Ding et al., 2020; Huang et al., 2021; Liu et al., 2021). The decrease of ESV we observed in the 
study area is similar to the findings found in other regions of Xinjiang. For example, ESV 
decreases have been reported in the Ugan-Kuga River Delta Oasis (Mayila et al., 2018), the 
Yarkant River Basin (Zhang et al., 2018), and the Ebinur Lake Nature Reserve (Zhang et al., 2019). 
However, ESV of the Yanqi Basin has increased evidently in recent years (Sidik et al., 2016). The 
urban population in the northern slope economic belt of the Tianshan Mountains is exceptionally 
dense, driving land use/cover pattern change and increasing pressures on environment. It is 
essential to reinforce the conservation of forestland, grassland, and waterbody in this region 
because of their relatively high ESV (Yuan et al., 2019). The ecological quality of these land 
use/cover types should also be monitored to protect ecosystems and enable them to sustain further 
population growth (Tang et al., 2017). 


5 Conclusions 


In this study, we used remote sensing images and several statistical data to analyze the 
characteristics of LUCC and ESV in the northern slope economic belt of the Tianshan Mountains. 
The composition of land use/cover types in the northern slope economic belt of the Tianshan 
Mountains changed significantly from 1975 to 2018. Land use/cover types showed different 
degrees of transformation over time, with the area of construction land showing the greatest 
increase. Specifically, the areas of cultivated land and construction land increased, while the areas 
of grassland and unused land gradually decreased. Although the areas covered by forests and 
waterbody increased from 1975 to 2005, they declined over the whole study period. The total ESV 
for the study region decreased by 4.00x 10° CNY from 1975 to 2018, mainly due to the changes of 
grassland ecosystems. The area of waterbody accounted for only about 3.50% of the total study 
area, but its CSS was high. Also, among the individual ecosystem service functions examined in 
the northern slope economic belt of the Tianshan Mountains, waste treatment, soil conservation, 
hydrological regulation, and biodiversity conservation were the most important ecosystem service 
functions. In conclusion, the results of our study can be used to promote sustainable development 
of the northern slope economic belt of the Tianshan Mountains, especially in integrating land 
management-ecosystem service approaches. 
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